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On quantum averaging, quantum KAM and 
quantum diffusion 

Sergei B. Kuksin 1 and Anatoly I. Neishtadt 2 

Dedicated to the memory of Mark Iosifovich Vishik. 
Abstract 

For nonautonomous Hamiltonian systems and their quantisations we 
discuss properties of the quantised systems, related to those of the corre- 
sponding classical systems, described by the KAM-related theories: the 
proper KAM, the averaging theory, the Nekhoroshev stability, and the 
diffusion. 

1 Introduction 

Consider a classical nonautonomous Hamiltonian system on the phase-space 
T*T d = R d x T d = {(p, q)} or T*R d = R d x R d with a Hamiltonian H(p, q, t): 

p = -V q H, q = V P H. (1) 

The corresponding quantum Hamiltonian operator is obtained by replacing in 
H(p, q,t) the variable qj, j = 1, . . . ,d, by the operator which acts on complex 
functions u(x) as multiplying by xj, and replacing eachp,,- by the operator f-jf-, 

where H is the Planck constant^ The Hamiltonian operator % = H(iV x ,x,t) 
defines a quantum system, and a classical problem of the quantum mechanics, 
streaming from its first years of existence, is to study (spectral) properties of 
the operator % and the properties of the corresponding evolutionary equation 

iKu(t,x) = Hu(t,x), (2) 

in their relation with the classical system (JTJ) . 
For example, if 

H(p,q,t) = \p\ 2 + V(t,q), (3) 
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then 



h 2 A + V(t, x) 



(4) 



i.e. % is the Schrodinger operator with the potential V. 

In this paper we discuss properties of the Hamiltonian operator Ji, corre- 
sponding to properties of system ([T]), described by the KAM- related theories. 
Namely, by the proper KAM, the averaging, the Nekhoroshev stability, and 
the diffusion (this list by no means is canonical; it corresponds to the authors' 
taste). We discuss results for quantum systems @ which we regard as parallel 
to the three classical theories above, mostly restricting ourselves to the case of 
periodic boundary conditions x £ T d and assuming that h = const. Scaling x 
and t in the dynamical equation ©, (Q} we achieve H = 1. A discussion con- 
cerning semiclassical limit h — > 0, when it is not appropriate to scale h to 1, 
is contained in Section [5] There we consider the equations in the whole space, 
x £ M. d , since for the periodic boundary conditions the corresponding results 
are less developed. 

All quantum results we discuss deal with non-autonomous equations @, 
((4]) , so their classical analogies are "KAM- related" theories for non-autonomous 
Hamiltonian systems ([3]). We do not touch very interesting, important and 
complicated problem of constructing eigenfunctions of nearly integrable Hamil- 
tonian operators by quantasing KAM-tori of the corresponding autonomous 
Hamiltonian systems (see [Laz93 ). 

Let u(t) be a solution of the equation ([2]), ((4]). Multiplying the equation 
by u and integrating over T d we get that |u(i)|| 2 = const. Write u(t,x) — 
^2 s u s (t)ip s (x), where {<p s } are eigenfunctions of the "unperturbed" Hamilto- 
nian operator. Then \ u s(t)\ 2 = const. What happens to the quantities |u s (t)| 2 
as t growths, i.e. how the total probability ^ |w s (£)| 2 is distributed between the 
states s £ 1 d when t is large? This is the question which is addressed by the 
theorems we discuss. 

Acknowledgments. The authors are thankful for discussions to Sergey Do- 
brokhotov, Hakan Eliasson, and Johannes Sjostrand. SK acknowledges the 
support of l'Agence Nationale de la Recherche through the grant ANR-10- 
BLAN 0102. 

2 Quantum averaging 

2.1 Averaging and adiabatic invariance 

Let a classical Hamiltonian ^ have the form 



where the unperturbed Hamiltonian \p\ 2 + V(t, q) , t — const, is integrable for 
each r. Let L, 1 < j < d, be the corresponding actions. The classical averaging 
principle (e.g., see in [AKN061 ILM88] ) implies that each action is an adiabatic 



H{p,q,et)=H e = \ P \ 2 + V{et,q) 



(5) 
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invariant, namely if u E (t) is a solution of the perturbed equation §\§h=h e , then 
Ij(u e (t)) stays almost constant on time-intervals of order e _1 . The averaging 
principle is a heuristic statement, and it does not always lead to correct results. 
The adiabatic invariance for classical systems is discussed in more details in 
Section [5] 

Now consider the quantum system, corresponding to the Hamiltonian above: 

it = Au + V(et,x)u), xeT 1 . (6) 

We assume that the function V(t,x) is C 2 -smooth bounded and denote by A t 
the linear operator in (jSJ), 

A t = -A + V(et,x). 

Let {(p s (r), s £ Z d } and {A s (r)} be the eigenvectors and the eigenvalues of A T , 
where each A s (r) is continuous in r. Let u(t, x) be a solution of ([5]), equal at 
t = to a pure state, 

u(0,aO = ¥>«,(()), (7) 
such that for each st, X So (et) is an isolated eigenvalue of A £t of a constant 
multiplicity. Consider expansion of u(t, x) over the basis {<Pb(t), s £ Z d }: 

u(t, x) = ^2 u s (t)cp s (et) . 

s 

The quantum adiabatic theorem says that u(t, x) stays close to the eigenspace, 
corresponding to A So (ei): 

Theorem 2.1. (M. Born, V. Fock WF28f and T. Kato IKat50f ) 

sup \ u s(t)\ 2 -> as e -> 0. (8) 

This is a very general result which remains true for systems in the whole 
space (when x € K d ) if the operator A et has mixed spectrum, but X So (et) always 
is an isolated eigenvalue of constant multiplicity, see in |LM88j . The case when 
this eigenvalues may be approached by other eigenvalues is considered in [AE99 . 

Both for classical and quantum systems, adiabatic theorems are considered 
often on infinite time interval — oo < t < oo under condition that the dependence 
of the potential V on time disappears fast enough as t — > ±oo, and the system 
is sufficiently smooth. In this case, for classical Hamiltonians with d = 1, 
difference of values of actions on trajectory at t — > ±oo tends to much faster 
than e as e — ► 0; in analytic case this difference is (9(exp(— const/e)), see |LL60j 
and references in |AKN06j . Sect. 6.4.5. For quantum systems, if for r — > — oo 
all the probability is concentrated in the states, corresponding to the eigenvalue 
X So (t), then all the probability but a very small remnant will be absorbed by 
these states as t — » +oo. In analytic case this remnant is 0(exp(— const/e)) 
|Nen93l IJP93] (this result also follows from the calculus, developed in |Sjo93| ). 

We will return to the quantum adiabaticity in Section [5] We note that also 
there are adiabatic theorems for systems where the Hamiltonian slowly depends 
not only on time, but also on a part of the space- variables, e.g. see [XKN06 , 
Sect. 6.4.1 for classical systems and [BDT06 for quantum systems. 



3 



2.2 Around Nekhoroshev's Theorem. 



Let us start with classical systems. Let H E (p,q) = h (p) + ehi(p,q), where the 
function h is analytic and steep (e.g., it is strictly convex, for the definition of 
steep functions see |Nek77j and [LM881 HHN06] ). Let (p(t),q(t)) be a solution 
of (JXJ) . Then there are a, b > such that 

\ P (t) - p(0)\ < e a V|i|<e £ ~\ (9) 

see in |Nek77, LM88, AKN06I. There are many related results. For example: 
let 

H e 0, q,t) = h lp) +eh 1 (cot-p,q), uj eM. N , 

where hi is an analytic function on T N x M d x T d , N > 1. Then for a typical 
uj estimate © is true. In particular, let us take 

H e (p,q,t) = \p\ 2 + eV{ujt;q). 

The corresponding quantised Hamiltonian is the operator —A + eV(u>t; x), and 
the evolutionary equation is 

u = Au + eV(u)t;x)u). (10) 

Do we have for solutions of (jTOJ) an analogy of the Nekhoroshev estimate ©? 
I.e. is it true that actions of the unperturbed system, evaluated along solutions 
of the perturbed equation (jlOp . do not change much during exponentially long 
time? It turns out that a weaker form of this assertion holds true, even when 
e = 1! Let us consider the equation 

u= Au + V(t,x)u), (11) 

and consider the squared r-th Sobolev norm of u: 

\\u\\ 2 r =Y^ k| 2 (l + | S | 2 r, rel. 

This is a linear combination of the actions for the unperturbed system with 
7 = 0. 

Theorem 2.2. (\Bou99dj ). Let V(t,x) = V(cjt,x), where uj e R N is a Dio- 
phantine vector and V is a smooth function on T N x T d . Then for each r > 1 
there exists c(r) such that any solution u(t) of (jTTJ) satisfies 

\\u(t)\\ r < Const • (lnt) c W||u Q || r , Vt > 2. (12) 

So if uq is smooth, then the high states u s stay almost non-excited for very 
long time. We miss a result which would imply that the quantity in calcu- 
lated for solutions of (|11[) . stays small for long time. 

It is surprising that a weaker version of this result holds for potentials V 
which are not time-quasiperiodic: 
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Theorem 2.3. ([Bou99b}) . Let V be smooth and C k -bounded uniformly in 
(t,x) for each k. Then for each r > 1 and a > there exists C a such that 

||u(t)|| r <C a t ||uo||„ Vt>2. 

Also see DellO . If the potential V(t,x) is analytic, then the norm ||u(i)|| r 
satisfies (fT!?)) . see fWan08]. We are not aware of any classical analogy of these 
results. 

3 Quantum KAM 

Let (p, q) £ M. d x T d . Consider integrable Hamiltonian ho(p) = \p\ 2 and its 
time-quasiperiodic perturbation H £ (p,q) = h (p) + eV(ujt, q), uj £ R n , where 
V is analytic. For the corresponding Hamiltonian equation we have a KAM 
result: For a typical (p(0),q(0)) and a typical uj the solution (p(t),q(t)) is time- 
quasiperiodic. 

The quantised Hamiltonian defines the dynamical equation (ITU1) . We regard 
the vector w as a parameter of the problem: cj 6 U (e K n . We abbreviate 
L 2 = L 2 (T d ,C) and provide this space with the exponential basis 

{e™- x ,s £ Z d }. 

For any linear operator B : L 2 — » L? let (B ao ,a,b € Z d ) be its matrix in this 
basis. 

The theorem below may be regarded as a quantum analogy of the KAM 
theorem above. For d = 1 it is proven in |BG01j . and for n > 2 - in FK09 . 
We do not know how to pass in this result to the semiclassical limit. 

Theorem 3.1. If e ^ 1, then for most u> we can find an ip-dependent complex- 
linear isomorphism ^((fi) — ^ e ^(tp), ip £ T N , 

$>(ip) : L 2 -> L 2 , u(x) i — y $>((p)u(x), 

and a bounded Hermitian operator Q = Q SiU , such that a curve u{t) £ L 2 solves 
eq. (1101) if and only if v{t) — ty(tu))u(t) satisfies 

v = i(Au — sQv) . 

The matrix (Q a b) is block-diagonal, i.e. Q ao = if \a\ ^ \b\, and it satisfies 
Q ab = {2ir)- n - d J J V{<p, x)e 1 ^ dxdtp + 0{e>), 7 > 0. 

Moreover, for any p £ N we have \\Q\\hp,hp < Ci and ||^(<^) —id \\hp.hp < eC^. 

Here "for most" means "for all u) £ U e C U, where mes (U \ U £ ) < e K for 
some k > 0" . In particular, for any w as in the theorem all solutions of eq. (| 10[) 
are almost-periodic functions of time. Their Sobolev norms are almost constant: 
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Corollary 3.2. For lu as in the theorem and for any p solutions of (JTUJ) satisfy 
(1 - Ce)|K0)||„ < \\u(t)\\ p < (1 + C7e)||u(0)|| p , Vt > 0. 
This property is called the dynamical localisaton. 

Proof. Since Q is block-diagonal, then ||v(i)|| p = const. Since v(t) = 
ty(t)u(t) and || if? — id \\hp,h p — £ C-2, then the estimate follows. □ 

Remarks. 1) Let n — 0. Then (|10p becomes the equation u — —i[Au + 
eV {x)uj . Theorem states that this equation may be reduced to a block-diagonal 
equation u — —iAu, where A a b = if \a\ ^ |6|. This is a well known fact. 

2) For n = 1 the theorem's assertion is the Floquet theorem for the time- 
periodic equation (|10p . In difference with the finite-dimensional case, this is a 
perturbative result, valid only for 'typical' frequencies w£l and small e. 

Proof of the Theorem. Eq. (|10[) is a non-autonomous linear Hamiltonian 
system in L 2 : 

S 11 

ii = -i—H e (u), H £ {u) = -(Vu,Vu) + -e(V{ip + tuj,x)u,u). 
ou 2 2 

Consider the extended phase-space L 2 x T™ x W 1 = {(u, tp, r)}. There the 
equation above can be written as the autonomous Hamiltonian system 

it = -i-^h £ (u,ip,r), 
ou 

if = V r h £ = lu, 
f = -Vtphs, 

where h e (u,ip 7 r 7 e) — to ■ r + i(Vu, Vu) + he(V((p, x)u, u). So h e is a small 
perturbation of the integrable quadratic Hamiltonian h = to ■ r + |(Vu, Vu). 
To perturbations of ho applies the KAM-thcorcm from EK10 . To show how 
this implies the Theorem 13. II let us write h £ as 

h e (u, f,r,e) = uj ■ r + — (Vu, Vtt) + ef(u, tp, r). 

In our case f — \ (V(<p, x)u, u) . The theorem below is the main result of |EK10j ). 

Theorem 3.3. There exist a domain O — {\\u\\ < 6} x T™ x {|r| < 6} and a 
symplectic transformation $ : O — > L 2 x T™ x W l which transforms h £ to 

h = uj 1 ■ r + - (Vu, Vu) + e{Qu,u) + f'(u, tp, r), 

where f = 0(\u\ 3 ) + 0(\r\ 2 ). 

Torus To = x T" x is invariant for the transformed system, so 'I* (To) is 
invariant for the original equation. This is the usual KAM statement. Now it is 
trivial since it simply states that u(t) = is a solution on the original equation. 
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But the KAM theorem above tells more. Simple analysis of the proof (see 
a Remark in [EK2]) shows that if the perturbation ef is quadratic in u and r- 
independent, then the KAM-transformations are linear in u and do not change 
uj. So the transformed Hamiltonians stay quadratic in u. Hence, the Hamilto- 
nian ho is such that /' = 0. That is, 

ho = uj' ■ r + ^{Vu, Vu) + e(Qu, u). 
This proves Theorem 13. II 

4 Quantum diffusion. 

Let (p,q) £ R d x T d . Consider H £ (p,q) = \p\ 2 + eV(ut, q), where uj £ R N and 
V is analytic. Then 

i) by KAM, for a typical uj and typical initial data (p , qo) the solution such 
that (p(0),q(0)) = (poiQo) is time-quasiperiodic; 

ii) for exceptional u> and (po,qo) we "should" have the Arnold diffusion: the 
action p(t) of a corresponding solution slowly "diffuses away" from po- 

As before, the quantised Hamiltonian defines the dynamical equation (flO]). 

Claim 4.1. Let d = 1, N > 2 and the potential V is nondcgenerate in a 
suitable sense. Then there exist a smooth function u(0,x) and ui £ R N such 
that 

limsup ||u(t)|| s = oo (13) 

t— >oc 

for some s > 1. 

An example of a time-periodic potential V, satisfying (|13|) . is given in |Bou99aj . 
It is conjectured by H. Eliasson that the validity of the Claim for a typical po- 
tential follows from the method of his work [Eli02 . Proof of this assertion is a 
work under preparation. 

5 Perturbed harmonic and anharmonic oscilla- 
tors. 

In Sections El 0] we deal with the evolutionary Schrodinger equation under pe- 
riodic boundary conditions. Some similar results are available for equations in 
the whole space with growing potentials: 

• Consider Schrodinger equation in Mr: 

ii = u xx + (x 2 + ^ix 2m )u + eV(toj, x)uj , 

where /i > 0, m £ N, m > 2; V(<p,x) is C 2 -smooth in tp, x and analytic 
in tp, bounded uniformly in tp, x. An analogy of Theorem 13. II holds. See 
|Kuk93j (Section 2.5) for the needed KAM-theorem. 
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• Due to Bambusi-Graffi BG01 , the result holds for non-integer m. That 
is, for equations 

u = — i( — u xx + Q(x)u + sV(ipo + tuj, x)u) , 

where Q(x) ~ |a;| Q ,a > 2 as \x\ — ¥ oo. The potential V may grow to 
infinity as \x\ — > oo. 

• Liu- Yuan [LY10J allow faster growth of V(x) in x. Their result applies to 
prove an analogy of Theorem 13. II for the quantum Duffing oscillator 

u = — u xx + x A u + exV((po + tuj, x)u) . 

• Due to Grebert and Thomann |GT11| . the assertion holds for the per- 
turbed harmonic oscillator 

it = — if — u xx + x 2 u + eV(ipo + tu>, x)u) . 
What happens in higher dimensions, al > 2 ? - This is completely unknown. 

6 Quantum adiabatic theorem in semiclassical 
limit 

In this Section we consider the classical system on T*M. d 
Hamiltonian 

H{p,q,r) = \p\ 2 + V(T,q), r = et, 
and the corresponding quantum system 

ihu = —h 2 Au + V(t,x)u = T-L T u. t — et, (15) 

(see ([4])). We assume that for each r the potential V(t, x) grows to infinity with 
|x|, so the operator T-L T has a discrete spectrum. 

We fix small enough e that allows to make some statements about the dy- 
namics of the classical system, and then pass to the limit as h — > 0. This limiting 
dynamics may be quite different from that in Section 12.11 when h is fixed and 
e — > 0, as it was demonstrated by M. Berry |Ber84j in the following striking 
example. Let d = 1 and potential V for r = const has two (non-symmetric) po- 
tential wells. Generically, for r = const and small enough h each well supports a 
family of pure quantum states localised mainly in this well. Consider a solution 
u(t, x) of equation (|15p with an initial condition which is a pure quantum state 
from the left well. For however small e there exists Hq = Ko(e) > such that 
if < h < ho, then for each t € [0, 1/e] the function u(t, ■) is localised in the 
same left well. On the other hand, under some rather general assumptions, for 
however small h there exist Sq — £o(h) and positive constants a\ < 02, such 
that ii < e < eo then the function u(t, ■) is localised in the right well for 
aih/e < t < ajije. 

Discussion of the case e ~ h is contained in |Kar90j . In what follows Eq, c, Cj 
are positive constants. 



= R d x R d with a 
(14) 
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6.1 Systems with one degree of freedom 

Assume first that classical Hamiltonian (fbi| has one degree of freedom. We sup- 
pose that V is C°°-smooth and that in the phase plane of the Hamiltonian sys- 
tem (H3J) for each r = const there is a domain filled by closed trajectories. In this 
domain we introduce action-angle variables I — I(p, q, r), x — xipi 1i T ) mod 2-7T 
(i.e. x G IT 1 )- Invert these relations: p — p(I, %, t), q = q(I, x, r). Suppose that 
there is an interval [ai,6i], < a\ < b±, such that the map I,Xi T ^ Pi<l, T is 
smooth for I £ [oi, &i], x S T 1 , r g [0, 1]. We express Hamiltonian (fH)) via the 
action variable and slow time: i?(p, q, r) = r). 

For e > let p(t),q(t) be a solution of the perturbed system with the Hamil- 
tonian H (p, q, et). 

Theorem 6.1. (see, e.g., \Am89p Th ere exist £o,c\ such that for < e < Eq 
we have 

\I(p(t), q{t),£t) - l(p(0), q(0), 0)| < eye for < t < 1/e . 

Now assume that for each r = const e [0, 1], and each /* € (oi, 6i) Hamil- 
tonian (fT4|) has a unique trajectory with the action 1 = 1*. Consider the 
corresponding quantum system (|15[) . The operator "H T has a series of eigen- 
functions ip s (r) = ip s (r, x) such that 

||¥>«(t)|| = 1, (p a (r,x)—t0 as x — > oo, (16) 

and the corresponding eigenvalues are A s (r) = E(I s ,t) + 0(h 2 ), where I s = 
h(s + l/2) e [oi, b\] (this is the Bohr-Sommerfeld quantisation rule, see |MF81j ). 
We assume that V is such that the convergence to zero in (|16p is faster than 
|x| in any negative power. Let u(t, x) be a solution of non-stationary equation 
(fT5|) with a pure state initial condition u{Q,x) — ip So {Q). Denote by P( a ,/3)(r) 
the orthogonal projector in L 2 (R) onto the linear span of vectors <p s { T ) with 
I s € (ct,f3). The approach in Bor78 leads to the following 

Conjecture 6.2. There exist So, c\ such that if < e < £o and < h < e, then 
for any m > 1 and a suitable Ca(m) > we have 

- ) • (17) 

Thus u(t, •) stays close to the eigenspace that corresponds to eigenvalues 
from O(e)-neighbourhood of E(I S0 ,et). 

6.2 Systems with several degrees of freedom 

Now let classical Hamiltonian (|14l) has d > 1 degrees of freedom. As before, 
we assume that V e C°°. For each r = const let the corresponding Hamil- 
tonian system be completely integrable and in its phase space there is a do- 
main filled by invariant tori. In this domain we introduce action-angle variables 
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I = I(p,q,r), x = x(PiQt t ) <= T d . Invert these relations: p — p(I,x,T~),q — 
q(I, x, t). Suppose that there is a compact domain A <e such that the map 
I,X, T ^ P:1t t is smooth for / £ A, x S T d ,r e [0,1]. We express Hamil- 
tonian (TH)) via the action variables and slow time, H(p,q,r) = E(I,r), and 
denote by uj(I,t) = dE/dl the frequency vector of the unperturbed motion. 
We assume that the system is non-degenerate or iso-energetically nondegen- 
erate (see definition in |Arn89] . Appendix 8). The dynamics of the variables 
(I,x){t) — {I,X){p{t)i<l(t), £ t) is described by a Hamiltonian of the form (see 
|Arn89j . Sect. 52F) 

H(I, X , r, e) = E(I, t) + eHi(I, x , r), (18) 

where Hi is a smooth function on A x T d x [0, 1]. 

Let Kq be a compact set in R 2d . For (po,qo) £ Kq denote by (p,q)(t) = 
(p,q)(t.,po-,qo) a solution of the perturbed system with initial condition (p,q)(0) = 
(Po,?o)- 

Theorem 6.3. (see, e.g., \AKN06l WM8^ ). If0<e<e , then 

I sup \I(p{i),q(t),et)-I(p(0),q(0),Q)\dp dqo<c ly fi. 

In systems with d > 1 degrees of freedom the value of action-vector as a 
function of time may change considerably for some initial conditions due to 
the effect of resonance between unperturbed frequencies, i.e. components of the 
vector uj(I, r). We say that there is a resonance for some (/, r) if (k-ui)(I, r) = 
for a suitable vector k £ 1 d \ {0} (here • denotes the Euclidian scalar product). 

Now consider corresponding quantum system (|15[) . Under some conditions, 
the operator W, T has a series of eigenfunctions <y5 s (r) = ip s (r,x), s £ Z d , satisfy- 
ing (fTo]l . with eigenvalues A m (r) = E{I mi r) + 0(h 2 ), where I m = h(m + jk) £ 
A, m £ Z+, and n £ Z d is the vector of the Maslov- Arnold indices |MF81j (the 
Bohr-Sommerfeld quantisation rule). Consider now the solution u(t,x) of non- 
stationary equation (|15j) with a pure state initial condition w(0,x) = tp mo (0). If 
we fix some small h and proceed to the limit as e — > 0, then Theorem 12.11 would 
apply. However, now we are interested in another limit, when a small e is fixed 
and H — > 0. Not much is known about the corresponding limiting dynamics. 
So we will formulate natural hypotheses about the limiting quantum dynamics 
as h — ¥ and will use them jointly with the known results about dynamics for 
classical Hamiltonian (|14[) with small e. 

For Theorem 12. ll to hold it is important that A mo (r) is an isolated eigenvalue 
for all r. Consider the distance between A m (r) and A TOo (r), where m,m £ Z d 
are such that m ^ tuq and |m — mo| ~ 1: 

A m (r) - A mo (r) = E(I m ,r) - E{I mo , T ) + 0(h 2 ) 

= {Im - I mo ) ■ w(/m„,T) + 0((I m - J mo ) 2 ) + 0{H 2 ) 

= h(m-m ) -w{I mo ,T) + 0{H 2 ). 



10 



Thus if there is no resonance at I mo , r, then distance between A mo (r) and nearby 
eigenvalues is ~ K. However, if there is a resonance k ■ w(/ mo ,r) = 0, then 
A mo +^fe( r ) — A mo (r) = 0(h 2 ) for integer v ~ 1. Thus classical resonances 
correspond to almost multiple points of the spectrum of the quantum problem. 
Therefore it seems that they should also manifest themselves in the quantum 
adiabaticity. 



For Hamiltonian (|14p there is a rather detailed information about dynamics 
in the two- frequency case d — 2. We will now use this information and the 
Bohr-Sommerfeld quantisation rule to state some conjectures about dynamics 
for the 2d quantum system (IT51) . 

Following P. Dirac |Dir25j we assume that 

dun dui 2 ( duji dui 2 \ dHi _ x 

W2 s7 " Wl a7 " {^W Wl aT J ^ > c (19) 

for all I, ip. General result by V. I. Arnold about averaging in two-frequency 
systems |Arn65, AKN06 implies that in this case 

\I(p(t),q(t),et)-I(p(0),q{0),0)\ < c ly fl for0<i<l/e. (20) 

On the basis of the Bohr-Sommerfeld quantisation rule and by analogy with 
Conjecture 16.21 it is natural to conjecture that for < t < 1/e the total prob- 
ability is mostly concentrated in the states, corresponding to actions 
from the Cy/e- vicinity of the original action I So . 

Now assume that instead of (fT9|) the following condition is satisfied (cf. the 
forth footnote) 

d^i duj 2 ^ _i , 01 v 
uj 2 — cji— — > c . (21) 

OT OT 

This is a particular case of a condition introduced by V. I. Arnold in |Arn65j . 
If, in addition to (|2Tj) . some general position condition is satisfied (see details 
in [AKN06 ), then estimate (f20| in which y/e is replaced with y/e\ lne| holds for 
all initial data outside a set of measure 0(y/e) |AKN06] . Sect. 6.1.8. The later 
set mainly consists of initial data for trajectories with capture into resonance, 
along these trajectories actions change by values ~ 1. Since for some initial 
data J(0),x(0) the solution I(t) is not localised in the vicinity of 1(0), then we 
should not expect for the quantum system (|15[) any estimate similar to that of 
Conjecture E21 where the amplitudes of eigenmodes tend to as h — > outside 
some small interval of actions. 

Consider classical Hamiltonian ((14)) under condition (|2T|) . Then the capture 
is only possible for a finite number of resonances, and the dynamics with a 
capture into resonance k\Wi + k 2 uj 2 = with co-prime k\,k 2 is the following 



4 Condition II19I I just means that the ratio of frequencies changes with non-zero rate 
along solutions of the system with Hamiltonian 1181 : u, |^(^) > c — 1 e. Similarly, con- 
dition 1211 1 means that ratio of frequencies changes with non-zero rate in adiabatic dynamics: 
^4( — ), . >c- 1 e. 

* at \ ujo ' I — const 
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|Nei05j . Denote (7,%)(t) = (I,x)(p(t),q(t),et). Suppose that at the initial 
moment t = we have no resonance: 

fciwi (7(0), 0) + fc 2 w 2 (J(0), 0)^0, 

and let t* G (0, 1) be the first moment when the resonance occurs: 

feiWi(/(0), r.) + fc 2 w 2 (I(0), r.) = 0. 
Then for < ei < r* the values of actions are approximately conserved: 

J(t) = 1(0) + 0(Vi hie). 

For t* < ei < 1 the system is captured into resonance, and evolution of actions 
is described by two relations: 

fciwi (J(i) , et) + fc 2 w 2 (/(*), et) = 0{ y/i In e), 
fca/^t) - kih(t) = fc 2 /i(0) - fci/ 2 (0) + 0(Viln £ ) . 

First of them means that the system stays near the resonance, while the sec- 
ond says that the dynamics has an approximate first integral. Jointly the two 
relations approximately define the trajectory I(t) for r* < et < 1. 

Based on this description and the Bohr-Sommerfeld quantisation rule, by 
analogy with Conjecture 16.21 we conjecture that for the quantum problem 
the capture in resonance in the classical system (I14[) results in transfer of an 
Ce-amount of the total probability from the vicinity of the initially excited pure 
state, corresponding to the action I so) to the vicinity of a state St £ Z 2 such 
that the lattice vector I{t) = h(s t + jk) satisfies the two relations above. This 
transfer happens for t > s~ 1 t*. When h— > 0, this Ce-amount stays positive of 
order e. 

For the dynamics of captured into resonances phase points also there is a 
more detailed description [Nei05 . Consider the resonant phase 7 = fciXi + ^2X 2 - 
It turns out that the behaviour of 7 is described by an auxiliary Hamiltonian 
system with one degree of freedom and the Hamiltonian of the form 

F = \fe (a(r)p2/ 2 + /( 7 , T ) + L(t) 7 ) . 

Here p 7 ,7 are canonically conjugate variables, function / is 2-7r-periodic in 7, 
and a, L ^ 0. In the phase portrait of the system for frozen r there are domains 
of oscillations of 7. Motion in these domains can be approximately represented 
as composition of motion along a trajectory of Hamiltonian F with frozen r and 
slow evolution of this trajectory due to the change of r. This evolution follows 
the adiabatic rule: the area surrounded by the trajectory remains constant. In 
the original variables p, q this motion is represented as a motion along slowly 
evolving torus. Angular variables on this torus are 7 and ip = l\(p\ + i 2 <^ 2 , 
where l\ and I2 are integers such that fc]7 2 — kiW = 1. This torus drifts along 
the resonant surface k\U\ + fc 2 w 2 = as it was described above. It is not known 
which quantum object corresponds to it. 
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